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On weak coupling expansion in models with unbroken symmetry * ' 
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An investigation of the weak coupling region of 2D SU(N) spin models is presented. An expansion of the free 
^\ ■ ener gy an d correlation functions at low temperatures is performed in the link formulation with periodic boundary 
conditions (BC). The resulting asymptotics is shown to be nonuniform in the volume for the free energy. 
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This paper deals with 2D SU(N) x SU(N) 
principal chiral model in the weak coupling re- 
gion whose partition function (PF) is given by 



J\DU X 



exp 



where U x € SU(N), DU X is the invariant mea- 
sure and we assume periodic BC. The main ques- 
tion discussed here is what is the correct asymp- 
totic expansion of nonabelian models at large (3 
and whether the conventional perturbation the- 
ory (PT) gives the correct expansion, the question 
addressed few years ago in |Q. As was rigorously 
proven, the conventional PT gives an asymptotic 
expansion which is uniform in the volume for the 
abelian XY model M. One of the basic theorem 
which underlies the proof states that the following 
inequality holds in the 3D XY model 



< exp(y/PM<l>x)) > < C, 



(2) 



where C is /3-independent and A(cj) + 2n) = (f>. <j) x 
is an angle parametrizing the action of the XY 



model S = J2a 



cos 



\- n ). It follows, at 



jx,n \ti Vx+r, 

large (3 the Gibbs measure is concentrated around 
<j) x giving a possibility to construct an ex- 
pansion around <p x — 0. This inequality is not 
true in 2D because of the Mermin- Wagner theo- 
rem (MWT) , however the authors of Q prove the 
same inequality for the link angle, i.e. 



< exp(yfj3A((/)i)) > < C , <f>i = <j) x - § x 



(3) 
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Thus, in 2D the Gibbs measure at large j3 is con- 
centrated around <pi ss and the asymptotics can 
be constructed expanding in powers of <f>i . In the 
abelian case such an expansion is equivalent to 
the expansion around <p x — because i) the ac- 
tion depends only on the difference <p x — 4>x+n arL d 
ii) the integration measure is flat, DU X = d(j) x . 

I'm not aware of any similar rigorous result 
proven for the nonabelian model. Nevertheless, 
because of the MWT one has to expect something 
similar to (||), namely 

< exp(V^argA(Tr^)) >< C , V t = U x U+ +n . (4) 

This is equivalent to the statement that the Gibbs 
measure at large (3 concentrates around Vj ~ /. 
That such (or similar) inequality holds in 2D 
models is intuitively clear and should also follow 
from chessboard estimates. It is assumed here 
that (U) is correct, hence Vi = I is the only sad- 
dle point for the invariant integrals. Thus, cor- 
rect asymptotic expansion should be given via ex- 
pansion around Vi = I, similarly to the abelian 
model. If the conventional PT gives correct ex- 
pansion, it must reproduce the series obtained by 
expansion around Vi = I. However, neither i) nor 
ii) holds in the nonabelian models, therefore it is 
far from obvious that two expansions indeed coin- 
cide. Our purpose here is to develop an expansion 
around Vi = I aiming to calculate the asymptotic 
series for nonabelian models. 

To accomplish this task we make a change of 
variables Vj = U x U+ +n in (0). PF becomes [| 



i exp 



fiJ2 Re TrVJ +ln J(U) 



(5) 
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j(v) = n 



l€p 



(6) 



Y\ p is a product over all plaquettes of 2D lattice, 
the sum over r is sum over all representations of 
SU(N), d r = Xr(I) and 

J] Vi = V n (x)V m (x + n)V+(x + m)V+(x). (7) 

JiV) is a product of SU(N) delta-functions 
which introduce the constraint Y\i ep Vi — I 3 . Its 
solution is a pure gauge Vi = U x ll x+n , thus two 
forms of the PF are equivalent. The correlation 
function T(x,y) =< Trf7 x Z7+ > becomes 



T(x,y) =<Tr ff Wi > , 



(8) 



where C xy is some path connecting points x and y 
and W; = Vi if along the path C xy the link / goes 
in the positive direction and Wi = VJ + , otherwise. 
The abelian analog of (0) reads 



n 



cos 4>i 



n e 



Jr(fi p 



(9) 
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</>p = (j) n {x) + (j) m {x + n)-(j) n {x + m)-4> m [x + n). 

The main step is to find a form for (|^) appro- 
priate for large-/? expansion. We consider SU (2) 
model and parametrize (a k are Pauli matrices) 



Vi = exp[icr fc w fe (/)] , ft = 1, 2, 3 , 



(10) 



1/2 



1/2 



,(11) 



where u> k (p) is a plaquette angle defined as 

V p = \{Vi=e^[ia k u k { P )]. (12) 



Then, the PF (§) can be exactly rewritten to the 
following form on a dual lattice (p — > ir) 



n 



2 L 



sin 2 W/ 



l e /3cosW, 



3 On the periodic lattice one has to constraint two holon- 
omy operators, i.e. closed paths winding around the whole 
lattice. Such global constraints do not influence thermo- 
dynamic limit in 2D. 



n 



w x 

sin W a 



n e 

x m{x) — — oo 



exp 



i / ctk(x)u)k(x) + 2Trim(x)a(x) 



k 



(13) 



where we have introduced the auxiliary field 
a k (x) and a(x) — Q2 k a k( x )) 1 ^ 2 ■ The deriva- 
tion of the PF ( p^|) is given in To perform 
large-/3 expansion we proceed in a standard way, 
i.e. first we make the substitution 

w fc (l) - (p)- x / 2 u k (l) , a k (x) -» (/3) 1 / 2 a fe (x)(14) 

and then expand in powers of fluctuations of link 
fields. Technical details are left to a forthcoming 
paper Q. The essential aspects are: 

1) . The main building blocks of the expansion 
are the following functions which we term "link" 
Green functions 

Gw — 2(5;. c — G x , x i — G x + n , x >+ n > + G x , x >+ n i 

+Gi+n,i', Dl(x ) = G XtX > — Gi+n.i'- (15) 

G x>x i is a "standard" Green function on the pe- 
riodic lattice. Unlike G XtX ', both Gw and D[(x') 
are infrared finite by construction. 

2) . Sums over representations are treated via 
the Poisson resummation formula. Due to this, 
zero modes decouple from the expansion after in- 
tegration over the auxiliary field. 

3) . An ensemble of fluctuations around Vi = I 
can be written as usual Gaussian ensemble ap- 
pearing due to integration over representations 
(over auxiliary field) in abelian (nonabelian) case. 

4) . Generating functionals (GF) appear to be: 



-hiGu'hp 



M (hi) = cxp 
for the XY model and 



1 



G u = 1 



(16) 



M(hi,s x ) = exp[ -s k (x)G x , x >s k (x') 
\a k {x)Di{x)h k {l) + \h k {l)G lv h k (l') } , 



(17) 



for the SU(2) model. h k (l) (s k (x)) is a source 
for the link (auxiliary) field and sum over all re- 
peating indices is understood. We list below some 
results proven in pj. 
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I. Asymptotics for the free energy of the XY 
model can be obtained from the expansion 



n 



00 1 d 2 



k=l 



M(hi) , (18) 



where Ak are known coefficients. It is easy to 
prove that this expansion coincides with the stan- 
dard PT answer. For example, the first order co- 
efficient is given by 



(19) 



II. Correlation function (g) in the SU(2) model 
up to the first order is given by 



T(x,y) = 1 



3 

4P 



E 



0(0- 



(20) 



where C xy is a path dual to the path C xy , i.e 
consisting of the dual links which are orthogonal 
to the original links I, I' <E C xy . The form of Gw 
guarantees independence of T(x,y) of the choice 
of the path C xy . After some algebra it is easy to 
recover the standard PT answer 

Y{x,y)=l-^D{x-y) , (21) 



L-l 



= 72 E 



1 — g 2 L'knX n 



fc„ = l 1J 2^in=\ Luh L Kn 



III. The main result of this study is the first 
order coefficient of the SU(2) free energy. The 
second order term in the auxiliary fields gives the 
following contribution at this order 



L 2 dlv ~ 8L 2 ^ * 

x.x' 



(22) 



Q XyX i is a second order polynomial in link Green 
functions (|l5|). Numerical study of ((2^) shows 
that j?C\ iv = b + ai \nL + O(^). All other 
contributions coming from the action, from the 
measure and from J(V) depend only on the link 
Green functions and arc finite. The final result is 



J2 C SU(2) = a o + ax In L + O(^) 
where at are (3- and L-independent. 



(23) 



Conclusions: The GF of the XY model de- 
pends only on Gw, guaranteeing the infrared 
finitness of the expansion. This is a direct conse- 
quence of the fact that the action of XY model is 
a function of <f>i. It is not the case in nonabelian 
models: the GF ([l7|) includes also dependence on 
G XjX ' creating potential danger for nonuniformity. 
The first order coefficient of T(x,y) is expressed 
only via Gw recovering thus the standard PT re- 
sult. A genuine nonabelian contribution appears 
only in the first order coefficient of the free enegy 
or in the second order coefficient of r(x, y). That 
there might be a problem with 2-nd order coeffi- 
cient of T(x, y) has been shown in jij, namely the 
PT with superinstanton BC (SIBC) gives a result 
different from the PT result obtained with peri- 
odic or Dirichlet BC. As was shown later, the PT 
with SIBC diverges at the third order ||. The 
interpretation has been given that it is just SIBC 
for which the asymptotics is nonuniform in the 
volume. Our result seems to support other in- 
terpretation: the correct asymptotic expansion 
of 2D nonabelian models is nonuniform in the 
volume. In particular, it means that the con- 
ventional PT does not take into account all the 
saddle points when the volume is increasing and 
becomes larger than the perturbative correlation 
length. We conclude thus either 1) to get the cor- 
rect expansion one has to expand a true function 
in the thermodynamic limit (like in ID) or 2) the 
inequality (Q) does not hold, e.g. the factor yffi 
is incorrect or 3) nonabelian models are not ex- 
pandable in 1//3 in two dimensions. 
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